On the universality of distribution of ranked cluster masses at critical percolation 
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The distribution of masses of clusters smaller than the infinite cluster is evaluated at the percola- 
tion threshold. The clusters are ranked according to their masses and the distribution P{M/L^ ,r) 
of the scaled masses M for any rank r shows a universal behaviour for different lattice sizes L {D is 
the fractal dimension). For different ranks however, there is a universal distribution function only 



in the large rank limit, i.e., P{M/L^ ,r) 



-vC 



g{Mr^ /L^) {y and (^ are defined in the text), where 



the universal scaling function g is found to be Gaussian in nature. 



Percolation is a classic example of systems with 
quenched disorder m. In a discrete lattice, sites or bonds 
are present with a certain probability and clusters are 
formed by connecting neighbouring occupied sites. At a 
critical probability, an "infinite" cluster appears for the 
first time which spans the whole lattice. 

The average mass or size of the spanning cluster is 
known to scale as M ~ L^ , where L is the lattice size 
and D the fractal dimension. In two recent papers [§,^, 
it was shown that when the clusters are ranked, the av- 
erage masses of the ranked clusters also show a similar 
scaling behaviour. This is true even for the clusters of 
large ranks, which are definitely smaller than the span- 
ning cluster. These clusters have been termed "effectively 
spanning" as their masses diverge with the lattice size 
although they do not really span the lattice. The be- 
haviour of the average scaled mass M/L^ as a function 
of the rank r was found to be 



(M/L^) - r 



„-A 



(1) 



where A can be expressed in terms of other known expo- 
nents of percolation as [|| 



A-l/(r-l). 



(2) 



Here r = 1 + d/D where d is the spatial dimension. It 
was also argued that the above behaviour is observed 
only in the asymptotic limit r — > oo. The (M/L^) vs. r 
curve actually changes its slope slowly (in a log-log plot). 
Hence for a given range of r, one can define an effective 
Ae//(r) with Ae//(r — > oo) given by (g). Very large rank 
would mean essentially clusters of size one or two in a 
finite lattice and these are not of present interest. 

Distribution function and its moments are useful for 
studying important properties of a system like multi- 
fractality, lacunarity etc. Distribution of the size of all 
the clusters, which is essentially the number of clusters 
of a given size (as a function of the size) in a dilute 
lattice, is well-known ||l| both at and away from criti- 
cality. Distributions of several other quantities like the 
size of the spanning clusters, chemical distances, short- 
est and longest paths on the percolation cluster etc. have 
also been studied in detail [Q-^. In general, at critical- 
ity, when the randomness is relevant, a universal non- 
Gaussian distribution function will exist Q. Although 



questions about distributions functions have been ad- 
dressed for quite some time, still a proper understand- 
ing is lacking in several areas [p],[7|-p|jll|] . Recently, the 
behaviour of the distribution of the largest cluster below 
criticality was also studied |12| . 

The existence of universal scaling functions for the dif- 
ferent quantities in the percolating lattice and the prop- 
erties of the ranked clusters inspired us to study the dis- 
tribution of the mass or size of these clusters at the perco- 
lation threshold. Although most of the quantities which 
have been studied earlier are directly related to the per- 
colating or spanning cluster (like the mass of the perco- 
lating cluster, the mass of the backbone, shortest path on 
the backbone etc.) the smaller clusters are no less impor- 
tant. In addition the remarkable fractal-like behaviour of 
the ranked clusters calls for further investigations. Our 
interest is particularly on the question of universality of 
the distribution function. 

In the simulation, the clusters in a square lattice 
(with helical boundary condition) are identified using the 
Hoshen-Kopelman algorithm. We rank the clusters at 
the percolation threshold irrespective of the fact whether 
the lattice is actually percolating or not. It may be noted 
that the ranked clusters may have degeneracy in the sense 
that there may be several clusters with the same rank in 
a particular realisation of the lattice. We checked, how- 
ever, that incorporating this degeneracy hardly affects 
the results. 

We first check the fact that the slope of the average 
cluster mass {M/L^) vs. r in a log log plot is indeed 
not unique in spite of M) and in agreement with |H]. We 
also verify that Xeff{r) has very weak finite size depen- 
dence, if any, as shown in Fig. 1. For r > 30, there is 
apparently some size dependence, but for the small lat- 
tices (e.g. L = 200,300 etc.), such ranks correspond to 
clusters which are not effectively spanning. Indeed, in 
H , the asymptotic value of A was found from very large 
lattices. However, one can obtain useful information as 
long as distribution functions are concerned even from 
relatively smaller lattices. 

The number of clusters of rank r with mass M/L^ 
is evaluated. The normalised probability distribution 
{P{M/L^ ,r)) of a cluster of scaled mass M/L^ and rank 
r is obtained by dividing this number by the total num- 
ber of clusters of rank r. This is shown for the ranks 
4,6,10 and 14 for several lattice sizes in Fig. 2. As in 



0, where distributions for the case r — 1 were consid- 
ered only, the bin sizes are proportional to 1/L , and one 
directly obtains a universal distribution for P{M/L^,r) 
for several values of L. Another interesting feature is, as 
one plots P{x — M/L^,r) for several ranks, it is found 
that the peaks of the distribution functions behave as 
Pmax{xmax) = x"^^ whcrc Xmax is the valuc of X at 
which the peak occurs. (This is shown by the straight 
line touching the peaks of the distribution in Fig. 2. in 
a log-log plot.) This behaviour of the peaks persist with 
a rank- independent value of C — 1-25 even for the higher 
ranks. 

We are more interested, however, in the behaviour of 
the probability distribution functions for different ranks 
for the same lattice size. The peak of the distribution 
P{M/L^, r) has a functional dependence on r as r^^''-'^^''^ 
from the above mentioned behaviour and eq p). How- 
ever, in general the behaviour of the entire distribution 
may not be as r'^^'ff^'^^ and we observe that it is better 
to expect a general form as 



P{M/L^,r)r~y'^ - g{Mry/L^) 



(3) 



when plotted against the natural scaling argument 
Mr^ /L^ . Here y is expected to be close to Xeff{r). The 
values of y are compared to different values of \eff{f) 
(corresponding to three different ranges of r) where the 
latter are obtained from a piecewise least square fitting 
of {M/L-^) vs r curves, g is a universal scaling function. 
Our attempt is to check whether one can actually obtain 
such a universal function for the distributions. 

While the data for the smaller ranges of r are taken 
from a system of L = 500, those in Fig. 5 correspond 
to that with L ~ 1000. The number of random config- 
urations generated are 10"* and 10"^ respectively for the 
two sizes. The results for the different ranges of r (as ap- 
propriate to the system sizes considered) are summarised 
below: 

Small r: For 4 < r < 14, we find that only one part 
of the curves (that beyond the peak value of P(a;, r)) are 
collapsing when plotted against the proposed scaling ar- 
gument with y ~ 1.75. Here the actual value of Ae//(r) 
is around 1.45. 

Intermediate r: For r values in a higher range (24 < 
r < 40), we find that the two parts of the curves collapse 
separately with different values of y; y ~ 1.4 for the por- 
tion beyond the peak, and y ~ 0.95 for the other portion. 
The value of Ae// for this range of r is found to be close 
to 1.25. 

Higher r: When one plots the scaled probabilities for 
even higher values of r, we find for the first time, a si- 
multaneous collapse of both sides of the curves with y 
between 1 and 1.1. The value of Ae// in this range is also 
~ 1.1 Hence a universal function is indeed obtained for 
large r values. We believe, for even higher ranges of r 
(for which reliable data can be obtained from larger lat- 
tices), the same behaviour will persist, with the value of 
y approaching the asymptotic value of A. Interestingly, 



for the smaller and higher ranges of r, y is neither equal 
to Ae// or the asymptotic value of A (at least for large x). 
However, in the scaling regime (i.e., for the large ranks), 

y - ^eff- ^ ^ 

The major portion of the universal distribution seems 
to fit well with a Gaussian distribution function of the 
form exp(— (x — xq)^ /a) with Q < x < oo, a ~ 0.0005, 
A -0.17 audio -0.1. 

Hence we obtain a distribution function in the follow- 
ing form 

P{M/L^,r) - r^^ exp{-{Mry/L^ - O.l)^)/^) (4) 

with y ~ 1.1 and C ci 1.25 for the higher ranks. 

Hence we obtain two most significant results in the 
present study: 

a) The exponent C — 1-25 for all ranges of the ranks. 
This is significant as while other properties of the system 
are rank dependent, this particular one remains constant. 

b) The existence of a Gaussian distribution: Most of 
the distribution functions studied earlier have yielded a 
more complicated universal function yj5|,|7|]. However, 
here also the data corresponding to very small values of 
Mr^ I L^ do not fall on the Gaussian fitting curve. 

It is difficult to relate Q to the known exponents in per- 
colation. Naively, if (Q) is to be derived from (Q), then 



(A//L^)=||lp(A//i^,r)d(-^) 



-A 



(5) 



gives C = 1.0 with y = X. This involves the approxi- 
mation that the mass of the cluster varies from zero to 
infinity. This approximation and also possible deviations 
from the Gaussian distribution may be responsible for the 
discrepancy between this value and the obtained value of 
(^; or it may simply be due to errors in numerical esti- 
mate. 

As already mentioned, the distribution for the prob- 
ability (per site) of clusters with s sites is known to be 
Q{s) ~ s""^ in a percolating lattice. One may expect 
that this behaviour can be extracted from P{M/L^,r) 
by calculating S^P(M/L^,r) as s = {M/L^)L^-'^, and 
a theoretical estimate of ( can be made. However, it has 
numerically been verified that one needs to include clus- 
ters of all ranks to obtain Q{s) in the above manner, and 
the absence of a universal scaling law for all r thus does 
not allow one to theoretically estimate C- 

In conclusion, the fact that the average cluster size ap- 
proached a rank independent scaling form given by (jj) 
only for large r, is consistent with our result that the 
universal form is obtained again only in the large r limit. 
One needs an exponent y to obtain a collapse of the data 
which should apparently equal A. However, y is greater 
than Ae// for the lower ranges of the rank. Surprisingly 
though, for the intermediate range of the rank, a data 
collapse is achieved for the smaller masses with a value 
of y very close to the asymptotic value of A. It is not 
clear how significant is this equivalence and whether it 
is purely accidental. In addition, we get an exponent ^ 



from the scaling behaviour of the probabiHty distribution 
which is independent of the rank. An approximate esti- 
mate of C is attempted to compare with the numerically 
obtained value. As in the cases of other quantities in 
percolation, here also a universal function is seen to ex- 
ist, which in contrast to the others is a simple Gaussian. 
The universal functions existing for each rank separately 
for several system sizes, however, have more complicated 
nature. 

The author is grateful to the computing centre of In- 
stitute of Physics, Bhubaneswar, where the programmes 
were run on a HPK9000/879 machine. She also thanks 
D. Stauffer for a critical reading of the manuscript and 
very useful comments. 



Fig 2: 2-cl, L = 200 to 500, distribution for tlie 4tli, Otii, lOtli and 14tli largest cluster 
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FIG. 2. The probability distribution for the ranked cluster 
masses are shown for ranks r — 4, 6, 10 and 14 (from right 
to left) for lattice sizes L = 200, 300, 400 and 500 against the 
scaled masses. The peaks of the distribution show a power 
law behaviour with M/L^ where the maxima occur for each 
rank. 
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4, 6, 10 and 14. The value of y is 1.75 here. 



Fig. 1 : average scaled mass vs. rank 
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FIG. 1. Variation of the average (M/L^) against the ranks 
are shown for lattice sizes L = 200, 300, 500 and 1000. Larger 
number of data (corresponding to higher ranks) are available 
for increasing lattice size. 





L = 200 O : 
L = .300 -1- : 
L = .500 □ . 
L = 1000 X 


- 


^*8 

^^0Xx - 
^Xx : 




■ 




0.01 0.1 

FIG. 4. The partial collapses of the data for the scaled 
distribution P{Mry/L'^) = P{M/L°)r-y'' for ranks 
r = 24, 28, 32 and 40 with two different values of y are shown 
separately. 
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FIG. 5. The collapse of the data for the scaled distribu- 
tion P{Mry/L'^) = P{M/L^,r)r-y'^ is shown for the ranks 
r = 70, 80, 90 and 100. The value of y is 1.1 here. 
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